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I This letter presents an adaptive spectrum sensing algorithm that detects wideband spectrum using 

sub-Nyquist sampling rates. By taking advantage of compressed sensing (CS), the proposed algorithm 

reconstructs the wideband spectrum from compressed samples. Furthermore, an £2 norm validation 

C/3 ' approach is proposed that enables cognitive radios (CRs) to automatically terminate the signal acquisition 

O 

once the current spectral recovery is satisfactory, leading to enhanced CR throughput. Numerical results 
show that the proposed algorithm can not only shorten the spectrum sensing interval, but also improve 

> 

. the throughput of wideband CRs. 
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I. Introduction 

Recently, cognitive radio (CR) has attracted much attention due to its capability of exploiting 
spectral holes and improving spectral utilization efficiency [[ll, This capability is fulfilled 
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by spectrum sensing which is defined as a technique for achieving awareness about the spectral 
usage and existence of primary users (PUs). With a "wider" spectral awareness, CR could exploit 
more spectral opportunities and achieve greater capacity. Thus, spectrum sensing over wideband 
spectrum becomes increasingly important for wideband CRs. 

To implement wideband spectrum sensing, CRs need some essential components, i.e., wide- 
band antenna, wideband radio frequency (RF) front-end, and high speed analog-to-digital con- 
verter (ADC). The wideband antenna and the wideband filter were well-developed as evidenced 
by [l3l and By contrast, the development of ADC technology is relatively behind: the 
achievable sampling rate of the state-of-the-art ADC is only 3.6 Gsps [l5l. To deal with this bot- 
tleneck, in the classic paper [6], Tian and Giannakis firstly applied compressed sensing (CS) [7J 
theory to CRs for acquiring wideband signals using sub-Nyquist sampling rates. Consequently, 
fewer compressed samples are required than predicted on the basis of Nyquist sampling theory. 
Furthermore, Wang et al. jHl proposed a two-step CS scheme for minimizing the sampling 
rate, where the actual sparsity was firstly estimated in the first time slot and the compressed 
measurements were then adjusted in the second slot. Additionally, Malioutov et al. @| studied 
a sequential CS approach where each compressed measurement was acquired in sequence. 

Against this background, the novel contribution of this letter is that an adaptive spectrum 
sensing algorithm is presented that utilizes CS theory to sense wideband spectrum by using 
an appropriate number of measurements. Different from the sparsity estimation scheme in [8], 
the proposed algorithm can adaptively adjust compressed measurements without any sparsity 
estimation efforts. Instead of the sequential measurement setup in [|9l, we acquire the wideband 
signals block-by-block from multiple mini-time slots, and gradually reconstruct the wideband 
spectrum using compressed samples until the spectral recovery is satisfactory. The remaining 
spectrum sensing time slots are utilized for data transmission, thereby enhancing the throughput 
of wideband CRs. Even with an unknown sparsity level, the proposed algorithm could still 
automatically terminate signal acquisition at the right time, leading to a robust spectral recovery 
as well as enhanced CR throughput. 

The rest of the letter is organized as follows. Section |ll] introduces the system model. Section 
Hn] proposes an adaptive spectrum sensing algorithm. Simulation results are presented in Section 
HVl with conclusions given in Section IVl 
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Fig. 1. Frequency and time frame in wideband CRs: (a) frequency frame, and (b) time frame. CR employs orthogonal 
frequency-division multiplexing techniques that divide the wideband spectrum into J subchannels. 



Suppose that CRs aim to exploit spectral holes within frequency band 0~ (Hz), as depicted 
by Fig. [Ha). Periodic spectrum sensing time frame is adopted as shown in Fig. [TJb) where 
~ r (second) is used for performing spectrum sensing and r ~ T (second) is reserved 
for transmitting data. During the spectrum sensing interval, all CRs keep quiet as enforced by 
protocols, e.g., at the media access control layer. Thus, the continuous signal received at the 
RF front-end of CR, i.e., Xc(t), is composed of only PUs' signals and background noise. By 
using sampling rate /a? over the observation time r, we could obtain a discrete time sequence 
x[n\ = Xc{j^), n = 0,1, ■ ■ ■ , A^— 1, in a vector form x E C^^^. Here, = r/jv is chosen to be 
a natural number. After spectrum sensing, CRs adopt orthogonal frequency-division multiplexing 
(OFDM) techniques that decompose the wideband spectrum into J orthogonal subchannels, each 
of which has bandwidth Bj = ^ (V j G [1, J]), as shown in Fig. [Ja). The subchannel index is 
denoted by j E [1, J] and PUs may present at any subchannels. For simplicity, let denote the 
set of subchannel indices where PUs present. However, based on the Nyquist sampling theory, the 
sampling rate is required to exceed 2W samples per second, i.e., /n > 2W; for a wideband CR, 
it leads to excessive memory requirement and prohibitive energy cost. This dilemma motivates 
us to employ CS technologies to reduce the sampling rate while retaining the spectrum sensing 
bandwidth W. 

CS theory indicates that, if a signal is sparse in some basis, it can be acquired by using a 
sub-Nyquist sampling rate; thus, fewer compressed samples are obtained than predicted using the 



II. System Model 
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Nyquist sampling theory. Mathematically, by using sub-Nyquist sampling rate fs (fs < 2Vr), 
the compressed samples y (y G C^^^^, M = rfs <^ N) can be written as 

y = (1) 

where $ denotes an M x measurement matrix. Notably, Tropp et al. [ fTO| cleverly implemented 
a CS system in which the measurement matrix is known and adjustable by changing pseudo- 
random sequences. For a comprehensive understanding of CS implementation, the reader is 
referred to ffTOl . 

In a CS-based spectrum sensing system, the goal is to reconstruct x or its discrete Fourier 
transform (DFT) spectrum X = Fx (F denotes a DFT matrix) from y. Then the traditional 
spectrum sensing algorithm, e.g., energy detection [fTTI . can be used to perform spectrum sensing 
using the reconstructed signal. For a robust signal recovery, the vector x is required to be sparse 
in some basis. Due to low spectral occupancy, it is believed that the received signal at CRs is 
sparse in the Fourier domain [6J. Thus, x is often assumed to /c-sparse (k < M <^ N) in the 
Fourier domain, which means that the DFT spectrum X consists of k significant components 
which are not negligible. If this spectral sparsity level k is known, we can choose the number 
of measurements M to secure the quality of spectral recovery, e.g., M = Cok\og{N/k) for a 
Gaussian measurement matrix, where Co denotes a constant [7]. Nevertheless, in a practical CR 
system, the spectral sparsity level is often unknown or difficult to estimate due to the dynamic 
activities of PUs. Furthermore, to avoid incorrect spectral recovery, traditional CS approaches 
tend to pessimistically choose Cq. Both phenomena can lead to more number of measurements 
and higher energy consumption, therefore, losing the advantage of using CS technologies. 

III. Adaptive Spectrum Sensing 
In this section, we study an adaptive spectrum sensing algorithm for wideband CRs. 

A. System Description 

Consider that a CS system, e.g., random demodulator [fTOl . is employed for implementing 
wideband signal acquisition as the discussions in Section |Ill Rather than sampling the wideband 
signal for the whole spectrum sensing interval in the traditional CS system, we propose to acquire 
the wideband signal step by step. The proposed algorithm aims to terminate the signal acquisition 
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TABLE I 

Adaptive Spectrum Sensing Algorithm. 



Initialize: Divide the spectrum sensing interval r into L mini time slots 
and set the mini time slot index I — 1 and accuracy e. 

Wliile the halting criterion is false and / < L, do 

a) . Acquire the compressed samples till the mini time slot /, 

resulting in the set of compressed samples yi. 

b) . Decompose the compressed samples yi into 

the training subset Ri and the testing subset Vi. 

c) . Estimate the wideband spectrum by applying a certain recovery 

algorithm to (|4j, leading to a spectral estimate Xi. 

d) . Calculate pi by using Vi and 

e) . If the halting criterion is true 

1) . Terminate the signal acquisition. 

2) . Perform spectrum sensing using Xi. 

3) . Choose frequency bands and start data transmission. 
Else: 1 = 1 + 1. 

Endlf 

Halting criterion: \pi/vi — 2S^\ < e. 

once the spectral recovery is satisfactory, and use the remaining spectrum sensing time interval 
for data transmission. The detailed algorithm is given in Table HI 

As shown in Fig. \V[h), the spectrum sensing interval is divided into L mini time slots where 
/ (/ G [1, L]) denotes the mini time slot index. Let yi (yi E C*^'^^) denote the set of compressed 
samples obtained from the beginning of spectrum sensing to the end of l-th mini time slot, and 
Ml denote the total number of measurements in yi, thus, < Mi < ■ ■ ■ < M^. Additionally, 
the sub-Nyquist sampling rate fs is chosen such that Ml = fs^ = C'o^maxlog(A^/fcmax) where 
^max denotes the maximum sparsity that can be estimated by long-term spectral observations. 
The set of compressed samples yi is then divided into two complementary subsets, i.e., the 
training subset Ri {Ri E C^^) for reconstructing the DFT spectrum, and the testing subset Vi 
(Vi E C^'^^) for validating the spectral recovery, where Mi = ri + vi. According to CS theory, 
the training subset and the testing subset can be written as 

Ri = ^ixi + n = ^iF^Xi + n (2) 
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and 

Vi = ^ixi+n = ^iF-^Xi + n (3) 

respectively, where F^^ is the inverse of DFT matrix, is a x measurement matrix, 
is a vi X N testing matrix, and n denotes the measurement noise modeled by circular complex 
additive white Gaussian noise (AWGN) with zero mean and variance 5^, i.e., n ~ CA/'(0,5^). 

By using a recovery algorithm, e.g.. Log-barrier approach [fT2l . we could obtain a spectral 
estimate Xi by solving the following problem: 

min llX^lli, s.t.: \\Ri - ^iF-'Xi^ < e (4) 

where e is a small recovery error threshold. Repeating this procedure, a sequence of spectral 
estimates, i.e., Xi, X2, ■ ■ ■ , Xi, will be obtained by increasing the total number of measurements 
Ml. Obviously, we would like to identify a "best" spectral estimate Xi that makes the spectral 
recovery error — X;||2 sufficiently small. If so, we can terminate the signal acquisition, and 
improve the throughput of CR system by using the remaining spectrum sensing time slots for 
transmitting data. However, the spectral recovery error \\Xi — Xi\\2 is typically unknown due to 
the unknown Xi when performing sub-Nyquist sampling. Hence, for a traditional CS system, 
the signal acquisition cannot be terminated at the right time. 

To identify the best spectral estimate, we propose to use the testing subset Vi for verifying 
the spectral estimate Xi. Specifically, we define the following verification parameter: 

pi = \\Vi-^iF-'Xi\\l (5) 

As we will see in the next section, if the verification parameter pi is close enough to 25'^ vi, the 
spectral estimate Xi is the best spectral estimate and the signal acquisition can be terminated. 

B. Performance Analysis and Comparison 

The termination metric in the preceding section is due to the fact that the best spectral estimate 
can be identified by validating the spectral estimate sequence and monitoring pi. 

Theorem 1: Let e > 0, g E (0, 1), and vi > 0. Given the sequence of spectral estimates 
Xi, - ■ ■ ,Xi, the best spectral estimate exists and is included in the sequence when the verification 
parameter pi satisfies 

PT[\pi/vi-26^\<e]>l- g (6) 
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where g = 2 exp 



2454+2(C/2+52)£ 



) 



, in which U denotes the measurement noise upper bound, 



i.e., n <U . 

The proof of Theorem 1 is given in the Appendix. 

Remark 1: It shows that, if the best spectral estimate exists within a given sequence of 
spectral estimates, the verification parameter should be within a certain small range around 
26^ with a high probability. This probability exponentially increases as the size of testing subset 
increases. In other words, if we monitor pi, we have a higher probability of identifying the best 
spectral estimate when using more measurements for validation. However, in another application 
scenario where the total number of measurements M/ is fixed, there exists a trade-off between 
training and testing. Even though allocating more measurements for validation (i.e., vi) achieves 
a higher probability of identifying the best spectral estimate, it leads to a degraded probability 
of successful spectral recovery because of fewer measurements for training (i.e., a larger vi leads 
to a smaller ri = Mi ~ vi). The investigation of this trade-off is an interesting issue for future 
research. 

Suppose that the signal acquisition is terminated at mini time slot /*, then the remaining 
time slots /* + 1, ■ ■ ■ , L could be used for transmitting data. Thus, the aggregate opportunistic 
throughput of the proposed CR system can be given by 



where Pf j is the probability of false alarm, Pj is the transmit power of CR transmitter, Hj denotes 
the magnitude channel gain between the CR transmitter and the CR receiver at subchannel j, 
and Nq denotes the noise spectral density. By contrast, a traditional CS system, the aggregate 
opportunistic throughput of CR system is given by 



It can be easily seen that the proposed system has superior performance than the traditional 
system due to /* < L in (|7]). 




(7) 




(8) 



IV. Numerical Results 



In simulations, we consider the following wideband signal: 




(9) 
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Fig. 2. Examples of: (a) the wideband spectmm X, and (b) the reconstructed spectrum X. The signal acquisition is terminated 
at mini time slot 10, with the number of compressed samples 2, 500. The SNRs of these 8 active subbands were set to random 
natural numbers between 7 dB and 27 dB. 



where sine fx) 



sin(7rx) 



a is a random time offset, z{t) is AWGN with zero mean and unit 



variance, and Ej is the received power at CR at subband j. The wideband signal consists of 
Nb = 8 non-overlapping subbands. The subband j is in the frequency range [fj — j-, /j + 
where the bandwidth 5^ = 10 ~ 30 MHz and the center frequency fj is randomly located in 

B B 1 

[-^, — ^] in which the overall bandwidth W = 2 GHz. The received signal-to-noise ratios 
(SNRs) of these 8 active subbands are random natural numbers between 7 dB and 27 dB. One 
time frame has length T = 10 /iS, in which the spectrum sensing interval is r = 5 fis. The 
spectrum sensing interval is divided into L = 20 mini time slots. Rather than using the Nyquist 
sampling rate /at = 2W = 4 GHz, we adopt the sub-Nyquist sampling rate fs = I GHz. 
The number of compressed samples in a traditional CS system is M = fsT = 5, 000, whereas 
N = /nt = 20, 000. The measurement matrix and the testing matrix follow the standard normal 
distribution with zero mean and unit variance. The measurement noise is assumed to be circular 
complex AWGN, i.e., n ~ CJ\f{0, 5^). The signal-to-measurement-noise ratio (SMNR) is set to 50 
dB. The energy detection approach in [[TT | is employed to detect PUs by using the reconstructed 
spectrum. For the data transmission, CRs adopt the transmit power Pj = 30 ~ 50 dBm. The 
channel between the CR transmitter and the CR receiver is assumed to be block slow fading 
channel with the path loss given by the 3GPP simulation guideline [fT3|: 127 + 301og^o(^)' 
where D (km) denotes the distance between the CR transmitter and the CR receiver. 

As we can see from Fig. |2l the wideband signal is composed of both high SNR subbands 
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shown. 

and low SNR subbands. Using the proposed algorithm, we can successfully reconstruct the 
wideband spectrum and terminate the signal acquisition at mini time slot 10, instead of mini 
time slot L = 20 when using the traditional CS algorithm. This is because, as shown in Fig. [3l 
the verification parameter becomes very close to 26"^ just when the (unknown) actual spectral 
recovery error becomes sufficiently small. Hence, if the result of Theorem 1 is used as the 
signal acquisition termination metric, the issue of excessive numbers of measurements can be 
solved. Fig. |4] shows that the wideband CR using the proposed algorithm outperforms the CR 
system using the traditional CS algorithms. The throughput gain improves as the transmit power 
increases. The reason is that, even with the same sub-Nyquist sampling rate, the proposed 
algorithm utilizes less time slots for performing spectrum sensing than that of traditional CS 
algorithms. 

V. Conclusions 

In this letter, we have proposed an adaptive spectrum sensing algorithm for improving the 
throughput of wideband CRs using CS technologies. It has been shown that the proposed 
algorithm can successfully reconstruct the wideband spectrum by using a few sub-Nyquist 
samples. Additionally, the wideband signal acquisition can be automatically terminated even 
if the actual spectral recovery error is unknown, thanks to the £2 norm validation approach. 
Furthermore, it has been proved that the proposed CR system can provide greater throughput 
than the CR system using traditional CS technologies. Simulation results have shown that the 
proposed algorithm can not only adaptively reconstruct the wideband spectrum by using an 
appropriate number of measurements, but also offer enhanced throughput for wideband CRs. 
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Appendix 
Proof of Theorem 1 

Due to the parameter setting in Section IIII-A[ the best spectral estimate (with zero or suf- 
ficiently small \\Xi — Xi\\2) should exist and be included in the sequence Xi,X2,-- - ,Xl. 
Assuming X;(/ G [1, -^^j) is the best spectral estimate, the verification parameter pi can be written 
as 



Pi 



\Vi-^iF-'Xi\\l = \\^iF-\Xi-Xi[ 



n\ 



\n\ 



(10) 



i=l 



where ^ and ri/ j denote the real and imaginary parts of the measurement noise, respectively. 
As riRi and n/_j are normally distributed with zero mean and variance 5^, we obtain E(n|.j) = 
E(n2 .) = and Var(n| J = ¥.{n\ --5'^f = Var(n| J = E{nli-6^f = 26\ Additionally, we 

Applying the Bernstein's inequality [fT4]| . we can obtain 



find Wn^i 



\i=l 

^ exp 
<1 exp 



=Pr 

e/2 



245%; + 2(f/2 + 52)^ 
Considering both (fTOl ) and (fTTI) . we obtain 



Pr [\pi-26\\ >l-2exp(- 



2A5\ + 2(t/2 + 52)^ 



(11) 



(12) 



Replacing by evi in ([T2)) . we complete the proof. 
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